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ABSTRACT 

Dynamic system models of linear and nonlinear hydrologic responses, 
as related to excess rainfall runoff processes, are employed and compared 
for flood hydrograph estimation, prediction and control. The linear 
response model is represented by the Dimensionless Instantaneous Unit 
Hydrograph and the nonlinear response model is based on a second order 
differential equation which includes the effects upon storage of unsteady 
flows and excess rainfall. Parameter estimations for these deterministic 
models are achieved with the application of direct search optimization 
techniques. A minimum variance Kalman filter algorithm is employed for 
the sequential estimation of the output and parameter of the linear 
model, in the presence of noise disturbances. In addition, a deterministic 
iterative estimator algorithm is used with the linear model in a recursive 
prediction of flows by way of minimizing the residual with respect to a 
dimensionless parameter as the control variable. The nonlinear model is 
examined on the basis of response surfaces for paired parameters and on 
the basis of the model output sensitivity to varying the parameters. 
Subsequently response noise models are determined and examined for 
possible typical real time model application. A satisfactory application 
of the nonlinear model was achieved in the prediction of the downstream 
hydrograph with input from excess rainfall and regulated dam releases. 
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FLOOD HYDROGRAPH PREDICTIONS AND CONTROL 

by , 

L.A. Logan 
INTRODUCTION 

Many of the hydrologic models that are currently available for 
flood flow predictions are mainly linear and deterministic with time 
invariant parameters. These models are often inadequate for some flow 
prediction and planning purposes because they do not, as a rule, include 
any explicit consideration of the uncertainties that are inherent in the 
form of a model and estimates of its parameters. In the real world, the 
rainfall- runoff process is truly nonlinear, nonstationary and is imbedded 
with both process and measurement noise disturbances. Although the 
probability characteristics of these noise processes may not be known or 
are known only partially, for practical purposes, they may be assumed to 
be zero mean Gaussian white noise processes (Jazwinski, (1970) ) . 

Early models, based on the Instantaneous Unit Hydrograph concept 
(lUH) (Nash (1957), Dooge (1959)) are linear and deterministic. A 
nonlinear lUH problem was treated by Singh (1964), and a reproduction of 
the excess rainfall-runoff process on the basis of a nonlinear reservoir 
storage model with time invariant parameters was achieved by Prasad 
(1967) . Amorocho (1967) , examined the general problem of nonlinear 
model predictions in the study of the rainfall-runoff process with 
special reference to the problems of predictive uncertainties. 

Recently, Hino (1973, 1975), Lee (1972), employed state estimation 
techniques of the modern Control Theory in the sequential estimation of 
parameters and states for linear and nonlinear hydrologic system response 
models. These approaches were attempts to account for some aspects of 
the stochastic nature of the response processes by including noise 
disturbances in solving the estimation problem. Duong, et al, (1975), 
and Szollosi-Nagy (1976), extended the Control Theory application in 
hydrology to include adaptive estimation and control of the input noise 
levels. 

This paper utilizes linear and nonlinear reservoir storage models 
in comparison for flood hydrograph predictions. Recursive filter 
algorithms are employed in the sequential estimation, with noise dynamics, 
of parameters and flows by the linear model. In addition, iterative 
estimation and prediction of flows are attempted by way of minimizing 
the mean squared residual with respect to the dimensionless model 
parameter as the control variable. The nonlinear model is tested for 
its applicability in the prediction of hydrographs with controlled input 
in the presence of excess rainfall. The objective is to implement 
simple but more realistic and effective models for flood hydrograph 
predictions and control. 
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Response Models 

The storage discharge relation for surface runoff processes, in the 
presence of disturbance effects on storage by unsteady flow and excess 
rainfall, is combined with the hydrologic continuity equation to produce 
a second order, nonlinear, differential equation representing a lumped 
nonlinear hydrologic response model: 

[a2D + N ai q^'-'-D + l] q = [l-a3D]u ... [l] 

Where q and u denote runoff output and excess rainfall input respectively; 

ai, a2 and a^ are storage coefficients for flow, unsteady flow and 

rainfall excess, respectively; 

N is a constant accounting for the nonlinearity of the system relative 

to lumped storage and discharge; and D E d/dt, is a differential operator, 

with dq/dt, defining the unsteady flow and du/dt, defining the rainfall 

intensity. 

By assuming linearity of the system (N=l) and negligible disturbance 
effect on storage, that is, uniform flow, dq/dt = 0, and uniform excess 
rainfall intensity, du/dt = 0, eq. [l] reduces to a linear, first order 
differential equation for the response model: 

[ajD + l]q^ = u^ ... [2] 

The solution of eq. [2] with unit step input, for equal n linear 
reservoirs in series and with the same time invariant parameter, results 
in the familiar Instantaneous Unit Hydrograph (lUH) (Nash 1957). The 
lUH may be reduced to the dimensionless lUH (Wu, (1963)): 



q 



Ft/tl n-l, -(n-l)(t/tp-l) ^^, ^3] 



where q„ and t„ denote peak flow and time-to-peak, respectively. Equation 
[3] respresents the unique functional for the linear model considered. 

A finite difference scheme may be employed in the numerical solution 
of eq. [1]. The resulting algorithm with rearrangement to solve for 
q(j+l) is: 



q(j + l) = 



u^^ o^^N-l^. 1-1 ^ Lj + D- ^(u(j + l)-u(j)) 
At2 2At J L ^^ 



- q(:)(l- ^^') -q(j-l)(!i-^q""'(j) 

At^ At^ 2At 

in which a central difference scheme is used for q, and a forward 
difference scheme for u. 



Dynamic System Models 

The equivalent dynamic system models may be represented by the 
following process and observation models (Jazwinski (1970)): 

dX(t) = f(X, U;t) + G(t) W(t), X{t^) = X^ ... [5] 
dt 

Z(t) = h{X,U;t) + V(t) ... [e] 

where X and U denote n-vector state output and p-vector input respectively; 

f (X,U;t) and G(t) denote n-vector functional characterizing the system 
and a noise coefficient matrix, respectively; W(t) is a p-vector process 
noise; 

Z(t) is a m-vector observed output; h(X,U;t) is a nonlinear functional 
characterizing the output and V(t) denotes a r-vector measurement noise. 

W(t) and V(t) are assvuned to be uncorrelated zero mean Gaussian white 
noise processes, with covariance matrices, Q and R, respectively. 

A recursive Kalman filter algorithm may be employed as an alternative 
technique to the conventional trial-and-error or least square methods in 
solving the estimation problem in real time application. For example, 
a discrete, minimum variance filter for the linear model, eq. L3j , 
based on eqs. [s] and [e] , comprises a conditional mean estimate (Jazwinski 
(1970), Sage and Melsa (1971)): 

A A '^ _ __ 

X(k|k) = X(k|k-1) + K[z(k) - HX(K|k-l)J ... [7j 
in which for the augmented system: 
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H =ri ^1 denotes the observation matrix; and K, denotes the steady 
state optimal Kalman gain matrix given by the Riccati equations: 

K = A P h"^ [h P H + r] ~ ... [9] 

and the steady state error covariance P, defined by: 

T T T r "I 

P = APA-KHPA + TQT . . . |_10j 

where A and Y denote the state coefficient matrix and noise coefficient 
matrix, respectively, for a discrete, linear, time invariant system. 

The sequential estimator results in A(k), H(k), F (k) and K(k) as 
updated matrices at each observation. Similarly, the covariance matrices 
Q(k) and R(k) , if required can be updated at each observation (Jazwinski 
(1970), Duong, et al (1975)). 



Noise Control Models 

In practice, the forms of the response models are usually inadequate 
and only estimates of the true values of model parameters are known; 
hence the residuals are often not truly white noise uncorrelated. A 
suitable noise model may be fitted to the residuals and used to modify 
or augment the response model in prediction. Improvement m predictions 
is achieved when the residuals of the augmented response noise model are 
approximately white noise uncorrelated. 

A suitable noise model of the residual may be incorporated into the 
model algorithm to provide optimal feedback control adjustments to the 
input ut as the control variable. The effect of the unobserved distur- 
bance, N^/ may be cancelled by setting (Box and Jenkins (1976)): 

Au^ = - v"-'"(B) (l-ct)iB)N^^ ... [ll] 

where Aut denotes the adjustment to be made on u^ at time t; and Nt+. 
denotes the model residual of the one-step-ahead prediction; V(B) is the 
impulse response function weights and <^iB denotes the lag-1 autoregressive 
noise model parameter. 

A suitable response noise model and feedback control algorithm may 
be developed for the linear or nonlinear model. 

APPLICATION AND DISCUSSION 

The response models are evaluated with experimental data in the 
Wilmot Creek IHD Representative drainage basin. Test applications were 
subsequently done in other selected project areas. Optimal estimates of 
the parameters for the deterministic forms of the models were determined 
with the use of direct search optimization techniques. "Golden ratio- 
technique was employed with the linear model and a "pattern search" 
technique was used with the nonlinear model. The minimum variance 
Kalman filter algorithm was also used in solving the linear estimation 
problem, with the consideration of associated process and measurement 
noise disturbances. Control noise model of the residuals was examined 
in relation to the nonlinear model. 

Linear Models 

As a result of the test applications in the Wilmot Creek basin, 
there appear to be sufficient variations in the model parameters, 
particularly parameter n, to indicate that these variations are due, 
probably, to parameter uncertainties as a result of influences by the 
process and measurement noise disturbances associated with the estimation 
problem. 

The recursive prediction with iterative estimation of the dimension- 
less parameter n as a control variable is demonstrated in Figure 1.1. 
The autocorrelation plots of the residuals indicated that the successive 
values of the residuals are approximately uncorrelated, in comparison 
with those of the predictions with the optimal value for n. The results 
of this control algorithms are satisfactory and although the method may 
be suitable for real time application, it is however, still dependent on 
the accuracy of the estimates of qp and tp to be determined by alternative 
empirical equations in relation to the accuracy of rainfall forecasts. 



Figures 1.2, 1.3 and 1.4 show the results of the application of the 
linear minimum variance filter for an alternative, one-step-ahead 
prediction with the consideration of input noise disturbances. The 
results of the sequential estimation of flows based on constant input 
noise levels are satisfactory. However, there is an appreciable improve- 
ment in predictions with an adaptive estimation of the input noise 
levels. The latter approach gave the more realistic results, in attempting 
to compensate for model error by way of adapting the noise levels to the 
residuals. In addition. Figures 1-3 and 1.4 demonstrate the stochastic 
nature of the sequential estimation as indicated by the propagation of 
the values for the steady state gains, K, and error covariances, P, 
respectively. In this example of the adaptive estimation, Q and R tend 
to converge very quickly to stable values from the initial guesses. 

The sequential minimum variance filter approach is to be favoured 
over the above deterministic recursive approach, seeing that in addition 
to providing simultaneous estimates of the state and parameter, it 
determines, also, estimates of the error covariances, from which reli- 
ability and accuracy of the predictions can be established. 

Nonlinear Model 

Preliminary analysis with the nonlinear model indicates that 
explicit consideration of nonuniform rainfall intensity, in solving the 
estimation problem, appears to be significant only in the larger catch- 
ments with a large a2 storage coefficient and in application with controlled 
inputs. Figure 2.1 shows examples of response surfaces for paired 
parameters. The parameters a^ vs a2 and N vs 0.2 display similar response 
surfaces with a distinct local minimum. However, parameters N vs a^, 
appear to have both a critical and sensitive relationship with respect 
to the objective function and describe a surface which appears to comprise 
contours approaching that of an elliptical hallow. 

The sensitivity of the model output to the varying of each of the 
model parameters, in turn, is demonstrated in Figure 2.2. The parameter 
02 influences and is directly proportional to the lag time to the peak 
flow. Both parameters aj and N influence the magnitude of the peak. 
For given optimal constants a2 and N, smaller values for ai are related 
to larger peaks and vice versa. However, for this example, in which the 
optimal ai<l and in relation to a small value for N, it is observed that 
with this combination for given constant 02 » a decrease in the value of 
N is associated with larger peaks and vice versa. However, in other 
applications, in which optimized parameter ai>l, a larger value of N is 
related to larger peak flows and vice versa. Figure 2.3(A) shows an 
example of the model application in a larger catchment, resulting in 
larger storage coefficients a^ and 02 • 

Figure 2.3(B) shows the result of the nonlinear model application 
in demonstrating its applicability in a typical problem solution. Here 
the upstream dam excess releases are regarded as controlled inputs and 
are used in addition to the excess rainfall and tributary inputs to 
generate the related downstream hydrograph. In this example the effect 
of storage coefficient, 03 , is included in the model calibration and 
prediction. Sensitivity analysis has indicated that a^ has direct 
influences on both the magnitude of the peak flow and time-to-peak 
{Figure 2.4(A) ) . 



Further applications, with the inclusion of the control noise model 

algorithm, results in significant improvement in the reproduction of the 

hydrograph, particularly in cases with less than satisfactory determin- 
istic calibration (Figure 2.4(B)). 

Table 1 summarizes the parameters and flow statistics for the 
nonlinear model application with a number of rainfall-runoff events. 

CONCLUSIONS 

The optimal values determined for the parameter, n, of the linear 
model, from its application to several storm events, exhibit significant 
variations even within a given catchment. Similarly, the hydrograph 
predictions with a constant optimal, n, gave less than satisfactory 
results. However, significant improvement in prediction is achieved 
with the use of a recursive and iterative estimation algorithm in which 
the model parameter, n, is treated as a control variable. Alternatively, 
the sequential estimation of the parameter, n, and prediction of flows, 
with the use of the minimum variance Kalman filter algorithm, attempted 
to account for imbedded noise disturbances in the solution of the problem. 

Further improvements in prediction are achieved as a result of 
attempts to compensate for some model error by the use of the residuals 
in an adaptive estimation of the noise covariances as part of the solution. 

The nonlinear response model applications resulted in satisfactory 
hydrograph predictions. However, the optimized parameters are, similarly, 
unique to a given storm event. The variations in the value of each 
parameter within a given catchment are probably due to the inherent 
system and observation noises associated with the modelling problem. 
There appears to be a consistent response surface relationship between 
paired optimized parameters, with the exception of paired parameters 
ai vs N. The optimization of the model is very sensitive to the proper 
combination of these parameters. The consideration of the effect of 
nonuniform input upon storage in the estimation problem appeared to be 
more significant in combination with controlled inputs. Satisfactory 
estimation of model parameters and hydrograph prediction are achieved 
with the combination of excess inputs from control dam discharges, 
tributary flows and excess rainfall. 

Examples of suitable response noise models are established for some 
selected catchment and storm situations. These could serve as preliminary 
models in simulation experiments for developing suitable control algorithms 
in attempts to compensate for model errors in the estimation problem. 

The sequential estimation and prediction algorithms are ideal for 
on-line implementation in solving real time flood forecasting problems. 
The combined nonlinear response and input control noise model, once 
suitably implemented for a particular problem site, could serve as a 
useful screening model in monitoring control of dam or reservoir releases 
in relation to any projected or target downstream rainfall-runoff hydrograph. 
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Table 1: Summary Statistics of the Nonlinear Model Applications. 



Basin 
(Area) 



Nithburg 
(326.3 km^) 

New Hamburg 
(551.7 km2) 



St. Jacobs 
(212.4 lan2)* 



West Montrose 
(194.3 }an2) 



Storm 
Date 



23/8/75 
#23/8/75 

23/8/75 
23/8/75 

16/5/74 

16/5/74 

#16/5/74 

16/5/74 

#16/5/74 



a 



10.2 
10.2 

13.8 
7.8 

6.6 
6.8 
6.8 

5.7 
5.7 



a- 



54.1 

59.1 

152.9 
145.6 

7.3 
8.0 
8.0 

0.41 
0.41 



-0.53 
0.83 

'V. 

-0.23 

-0.23 

1.0 
1.0 



N 



0.7 
0.7 

0.05 
0.10 

0.56 
0.55 
0.55 

0.48 
0.48 



±L 



-0.03 



0.83 



0.30 



Se 
(mm) 



0.13 
0.12 

0.15 
0.14 

0.64 

0.63 
0.39 

0.34 
0.31 



qp (mm/h) 



me as. 



2.3 
2.3 

1.4 
1.4 

9.0 
9.0 
9.0 

11.7 
11.7 



cal. 



2.1 
2.1 

1.2 

1.2 

6.6 
6.6 
7.8 

11.0 
11.1 



t (h) 
meas . cal 



16 
16 

22 
22 

19 
19 
19 

18 

18 



16 
16 

23 
23 

20 
20 
19 

19 
18 



Depth (mm) 
meas . cal , 



40.0 
40.0 

27.7 
27.7 



146.3 
146.3 



39.2 
39.4 

29.0 
29.0 



91.9 86.1 
91.9 85.6 
91.9 93.8 



144.7 
144.7 



ai - storage coefficient - flow 

02 - " '* - unsteady flow 

ag - " " - input 

N - nonlinear parameter 

^l - noise-model pareimeter 

S - standard error of estimate 



Nith/Nithburg (02GA038) 
Nith/New Hamburg (02GA018) 
Conestoga/St. Jacobs (02GA006) 
Grand/West Montrose (02GA034) 



S 



- peak flow 



tp - time-to-peak flow 

vol - volume of surface runoff 

'^ - assumed zero for 03 

* - drainage area: partial, up to the dam location 

# - with control noise-model 



Southern Ontario Streamflow 
gauging stations 
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Figure 1.1: Hydrograph Prediction: Linear Model (DIUH) 
(Wilmot Creek Basin - Event 16/8/69). 
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Figure 1.2: Hydrograph Prediction - Sequential Filter Estimation: 

Linear Model (DIUH) (Wilmot Creek Basin - Event 16/8/69) 
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Propagation of Error-Covariance, Pji* for Sequential 
Filter Estimation of Flow, q; Linear Model (DIUH) 
(Wilmot Creek Basin - Event 16/8/69). 
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Figure 1.4: Propagation of Kalman Gain, Ki, f or Sequential Filter 
Estimation of Flow, q; Linear Model (DIUH) (Wilmot 
Creek Basin - Event 16/8/69). 
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Figure 2.1: Parameter Response Surfaces: Nonlinear Model (Hflmot Creek 
Basin - Event 20/9/71). 
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Flp.ure 2.2: Kydrograoh Sensitivity to Varying Pariinieturs: Nonlinear Model (Uilmot Creek Basin - Event 21/9/71). 
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Ftgure 2.3: Hydrograph Predictions: Nonlinear Model Application (Grand River Basin 
23/8/75. (B) West Montrose 16/5/74). 
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Figure 2 . A : HyJrograph Prediction: Nonlinear Model Application - (A) Sensitivity to Varying Parameter 
03 - (B) Effect of Input Control Noise-Model. 
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